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ABSTRACT 

During a Type-I burst, the turbulent deflagation front may excite waves in 
the neutron star ocean and upper atmosphere with frequencies, ~ 1 Hz. These 
waves may be observed as highly coherent flux oscillations during the burst. 
The frequencies of these waves changes as the upper layers of the neutron star 
cool which accounts for the small variation in the observed QPO frequencies. In 
principle several modes could be excited but the fundamental buoyant r— mode 
exhibits significantly larger variability for a given excitation than all of the other 
modes. An analysis of modes in the burning layers themselves and the underlying 
ocean shows that it is unlikely these modes can account for the observed burst 
oscillations. On the other hand, photospheric modes which reside in a cooler 
portion of the neutron star atmosphere may provide an excellent explanation for 
the observed oscillations. 

Subject headings: stars : neutron, oscillations - X-rays : bursts, binaries 



1. Introduction 

As material accretes onto the surface of a star, the generation of nuclear energy may 
be stable or unstable depending on the rate of accetion and the properties of the underlying 
star. Unstable nuclear burning on the surface of a neutron star manifests itself as Type-I 
X-ray bursts (Hansen & van Horn 1975; Grindlay et al. 1976; Woosley & Taam 1976; Joss 
1977; Lamb & Lamb 1978; Lewin et al. 1993; Strohmayer & Bildsten 2003, the final two are 
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reviews) . The accumulation of material on the surface of the star is hkely to be sufficiently 
chaotic that the nuclear burning ignites at a particular point, so one would expect that as 
the burning envelopes the stellar surface that the observed flux to vary at approximately 
the spin frequency of the star (Joss 1978). Spitkovsky et al. (2002) and Zingale et al. (2001) 
present recent models for the growth of the flame other the stellar surface. The discovery of 
these oscillations had to wait for the launch of RXTE (Strohmayer et al. 1997). Surprisingly 
the frequency of the observed oscillations varied during the burst which Strohmayer et al. 
(1997a) argue is a hallmark of the conservation of angular momentum during the expansion 
and contraction of the atmosphere. 

The detailed models of Gumming & Bildsten (2000) estimate the radius expansion 
expected during the burst. Heyl (2000) calculated the general relativistic corrections required 
to translate the predictions of Gumming & Bildsten (2000) into observable quantities, and 
found that the observed frequency shift would be 30% - 50% of that predicted by Gumming 
& Bildsten (2000). Although Abramowicz et al. (2001) (and later Gumming et al. 2002) 
found a error in the derivation by Heyl (2000), Abramowicz and collcages also found that 
the general relativistic decrement is large. Gumming et al. (2002) also discovered an error in 
their earlier work (Gumming & Bildsten 2000) which lead them to overestimate the frequency 
shift by a factor of two even in the Newtonian context. 

Gumming et al. (2002) pointed out that the frequency shifts observed by Wijnands 
et al. (2001) and Galloway et al. (2001) were too large to be accounted for by the radius 
expansion models of Heyl (2000) and Gumming & Bildsten (2000). This letter examines an 
alternate model for both the observed frequencies during Type-I bursts and their evolution 
(first suggested by Gumming & Bildsten 2000). The burst may excite modes in the neutron 
star ocean or photosphere that exhibit themselves as dark and light regions on the surface. 
Because these regions are associated with ocean and atmospheric waves, they move relative 
to the surface of the star. The properties of the atmosphere change as the burst subsides, 
so the natural frequency of the modes shifts accounting for the observed frequency shifts. 

The following section of this letter identifies the modes excited by the burst and esti- 
mates the observational footprint of these modes. The last section examines the implications 
of these results. 



2. Waves 

Neutron stars exhibit several modes of oscillation (McDermott & Taam 1987; Lee & 
Strohmayer 1996). The modes that may be excited by the motion of the deflagation front 
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across the surface of neutron star undergoing a Type-I X-ray burst and may explain the 
observations of the burst oscillations share the following features. 

• They have a period on the order of one second. The observations find that the 
observed frequency shift is a few cycles per second. A mode frequency {u) of several 
Hertz is well matched to the timescale for the defiagation of the accumulated fuel to 
envelope the stars ~ 1 s (Strohmayer et al. 1997b), so it is unlikely that much slower 
modes would be excited. The spin frequencies of the stars {il) are several hundred 
Hertz so g = 2Q/uj ~ 10^ 

• The modes should travel westward. Here, m > denotes a westbound mode. The 
observed frequency is assumed to be slightly less than the spin frequency of the star, 
and it increases as the observed photon spectrum cools and dims; consequently, the 
modes should travel westward, i.e. in the opposite sense of the star's rotation. Some 
observations have found the opposite trend during a portion of the burst, so eastbound 
modes may also be interesting. 

• The mode should have no latitudinal nodes, or = 0. The modes of a rotating 
star are squeezed near the equator. Except for special geometries, a band both above 
and below the equator is visible throughout the star's rotation, so if the mode has 
latitudinal modes, much of the variability will be averaged out. 

• The azimuthal eigenvalue (m) should be 1 or -1. Observations have generally 
found that observed oscillation has a frequency of approximately the spin frequency 
of the star. Miller (1999) found that 4U 1636-536 may be exceptional, so the |m| = 2 
case will also be discussed. 

• The main observed mode should have no radial nodes in the radial displace- 
ment, or n = 1 {n counts the number of radial nodes in the transverse displacement). 
This ensures that any modes with a similar angular dependence will be well sepa- 
rated in frequency, supporting the argument that only a particular mode is excited 
and observed. 

Three prime candidates are the g^— modes, the buoyant r— modes and the Kelvin modes, 
of the neutron star ocean and upper atmosphere. In a rotating neutron star, the modes 
with various values of numbers of azimuthal nodes (m) are not degenerate in frequency. The 
frequency of the modes depends strongly on whether the electrons contribute to the entropy 
of the gas. If they are degenerate, the Brunt- Vaisala frequency (A^) is reduced by the ratio of 
the ion pressure to the total pressure; however, this reduction is mitigated by a corresponding 
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increase in the pressure-scale height (H). A first approximation to the frequency of a surface 
wave is 

/---. (1) 
27rR ^ ^ 

I shall examine modes in a semi-degenerate material with varying contributions of ion pres- 
sure to the total pressure a = Pi /Ptotai and non-degenerate material with varying contri- 
butions of gas and radiation pressure, /3 = -PgasZ-Ptotai- The heat transfer is assumed to be 
dominated by photons scattered by electrons in the nondegenerate regime and by electrons 
scattered by ions in the degenerate regime. 

The results of Bildsten et al. (1996) provide a touchstone for the various results. Extend- 
ing their results to a scmidegenerate regime yields the following estimates for the frequency 
of the modes in the rotating frame of the surface of the star, 



/...o - ^ = 2.37Hz 2A^^^^-^ — 
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Here A is the mean atomic weight of non-degenerate species (here, the nuclei in the ocean), 
Pb and pt are the densities at the top and bottom of the excited layer (Bildsten & Cutler 
1995, give the dependence on ocean depth), and n > is the number of radial nodes. If one 
takes the limit as a vanishes, the Bildsten et al. (1996) result obtains, and as a approaches 
unity for a fixed value of A the frequency estimate actually decreases slightly (by about 
40%). The frequency of the g-mode is proportional to the Brunt 

The ocean lies below the burning layers from p ~ lO''' g cm~^ to p ~ 10^ g cm^^, 
yielding 4.2 for the coefficient of is Eq. 2. The frequency of the first harnomic is about 
half that of the fundamental. The depth of the ocean depends sensitively on the charge of 
the nuclei produced by the burning. The unstable burning during Type-I bursts produces 
mainly iron group elements (Lewin et al. 1993), so the ocean is dramatically more shallow 
than in the higher accretion rate Z sources where the ocean consists mainly of CNO elements 
and extends to p ~ 10^^ g cm-=^ (Bildsten & Cutler 1995). 

In a nondegerate layer, the mode frequencies are several times larger. 
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where fi is the mean molecular weight of the material (this includes the electrons because in 
the non-degenerate case they contribute to both the entropy and the pressure) and /3 is the 
ratio of the gas pressure to the total pressure. To obtain these expressions, thermal bouyancy 
is assumed to dominate (Cumming & Bildsten 2000) and the expressions of Bildsten & 
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Cutler (1995) and Bildsten et al. (1996) have been scaled by the ratio of the Brunt- Vaisala 
frequency in the nondcngcrate atmosphere (Bildsten 1998) to the degenerate ocean. If one 
takes j3 ^ 1 and a — > 1 in expressions Eq. 2 and 3 for the same composition, the two 
expressions agree. The increase in the frequency estimates results entirely from the decrease 
in the mean molecular weight of the nondegenerate species, because the electrons are no 
longer degenerate. 

Both equations 2 and 3 include a logarithmic factor involving the density of the top 
and bottom of the layers. Bildsten et al. (1996) derived this factor when they considered a 
degenerate layer. It is unclear whether it is indeed appropriate in a non-degcncratc layer. 
This factor changes the estimate of the mode frequencies by at most a factor of two for the 
situations considered here. Furthermore, in the non-degenerate regime the mode period is 
shorter than the thermal time for the layer, so one would expect the modes to be isothermal 
rather than adiabatic. Bildsten & Cutler (1995) examine this point in their appendix and 
find that the lack of adiabaticity changes the structure of the mode and most importantly 
the assumed energy in a photospheric mode. However, they did not determine how it would 
change the frequency of a mode restricted to an isothermal region. 

Unless the temperature is small or n is large, the frequency of the oscillation will be 
much larger than one Hertz. The latter possibility is disfavored according to the criteria 
listed above. Furthermore, to meet the requirements above, A ~ 1. 

Longuet-Higgins (1968) derives the mode structure for the ocean of a rotating fluid shell 
which includes some modes not discussed by Bildsten et al. (1996). As q increases, for the 
gi— modes, A approaches 

A~(2i/ + l)V (m > 0,1/ = 0,1,2,...). (4) 

where z/ > 1 is the number of latitudinal nodes in the pertubation to the northward velocity 
of the fluid. 1^ = u + 1 for y-modes. 

Bildsten et al. (1996) also discuss a set of eastbound modes which Longuet-Higgins 
(1968) identifles as the Kelvin modes with 

A ~ (m < 0, z/ = 1) (5) 

which have l/^ — 0. 

For the buoyant r— modes (which all travel westward), A approaches, 

A~m2(2i/ + l)-2 (m> 0,1/ = 1,2,...). (6) 

and L, — — 1 for the r-modes. 
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As e = (fx increases the modes become more and more localized near the stellar equator, 
and the eigenfunctions resemble parabolic cylinder functions which may be expressed as the 
product of a Hermite polynomial and a Gaussian. The vertical displacement (or equivalently 
the pressure pertubation) for the mode ^ is given by the following expressions as e goes to 
infinity (Longuet-Higgins 1968), 

CKelvin ~ l,l/4e-|'^%-^+-* (8) 

m 

where t] = e^/^ji and ji = cos 9. Most of the excitation in a mode lies between < r] < \/v. 
For the different modes, this leads to the following range in /x, 

{q~^ (2 + v~^Y^^'^ gi-modes 
g-i/2|^|-i/2 Kelvin modes. (10) 

^-1/2 1 I -1/2 (^^2 _|_ ^-^1/2 ]-,yQyg^n^ r- modes 

The buoyant r-modes occupy the widest band near the equator; the variability is approx- 
imately proportional to the square root of the width of the band, so the r— modes should 
exhibit the largest variability. 

The next subsection focuses on the modes which are likely to be excited in the ocean and 
what their observational footprint would be. The modes tend to become more localized in the 
equatorial regions as the spin of the star increases (Bildsten et al. 1996). The observations 
indicate that the excited modes typically have |m| = 1 (Strohmayer et al. 1997a) or possibly 
|m| = 2 (Miller 1999). 

The frequency of the modes in the observer's frame is given by a;/ = a; — mil (neglecting 
the gravitational redshift for now), so for m > one would see the frequency of the vari- 
ation increase as the temperature of the ocean decreases. This is what is usually observed 
(Gumming & Bildsten 2000); therefore, the detailed calculations will focus on the m — \ 
and m — 2 modes with frequencies of about one cycle per second. 

2.1. Visibility 

Galculating the bolometric pulsed fraction for the mode provides an estimate the vis- 
ibility of the various modes. Because the frequency of the mode is about one Hertz and 
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the rotational frequency of the star is about 300 Hz, a value oi q = 300 is appropriate to 
calculate the modes. The variability of these modes decreases as q~^^^. The observer lies at 
a latitude of 30° and compare the total flux observed when the observer lies above a bright 
patch (intensity maximum) patch to when the observer lies above a faint patch (where the 
intensity vanishes). The pulsed fraction so defined is 

_ /(bright) - /(faint) 
/(bright) + /(faint) 

and can range from —1 to 1. The pulsed fraction is negative if the total flux observed from 
the star is larger when the observer lies above a faint region. 

Lacking a specific model for how the presense of a mode induces variations in the fiux, the 
effective temperature of each surface element is assumed to vary as the linear combination of 
a constant and the vertical displacement induced by the mode of interest. The normalization 
of the mode is selected so that the effective temperature varies over the surface from zero to 
twice the underlying value. Because the observed variability is small during the tail of the 
burst, the particular choice of coupling is unimportant; one can always choose a different type 
of connection between the mode and the observed flux and rescale the size of the underlying 
wave. 

The gravitiational defocussing of the neutron star surface in the Schwarschild geometry 
{e.g. Heyl & Hernquist 1998; Page 1995) and limb darkening also affect the variability. The 
limb darkening is modeled by allowing the intensity from a given surface element to vary 
as the cosine of the zenith angle {e.g. Perna et al. 2001). Figure 1 presents the results for 
the first few modes with m = 1 and m = 2. The inclusion of limb darkening increases the 
expected variabihty significantly and also alleviates the effects of gravitational defocussing. 
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Fig. 1. — Pulsed Fraction for various modes as a function of M/R for an observer latitude 
of 30° (left) and 60° (right). The upper curves include a rudmimentary treatment of limb 
darkening. A negative pulsed fraction indicates that the star is brightest when the observer 
is directly above a dark region. R^^ denotes the Rossby mode with the particular value of 



u and m. denotes a Kelvin mode. 
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As expected from the earlier discussion, the buoyant r— modes exhibit the largest vari- 
ability. Modes with a larger number of nodes v show less variability as do modes with 
\m\ > 1. The Kelvin modes are more concentrated near the equator, and therefore exhibit 
less variability. The (7— modes which arc extremely concentrated near the equator show neg- 
ligible variability when averaged over the visible portion of the star with a static background. 
Of all the modes in the ocean for a given excitation, the lowest-order r— mode exhibits the 
largest variabihty and also satisfies the observational constraints. The absolute variability 
of a particular mode depends on the local coupling with the mode and the emergent flux; 
the details of this coupling are beyond the scope of this paper, but the relative variability 
between the different modes will not depend on these details as long as the couphng is local. 



2.2. Veir lability Frequency 

With a mode excited, a distant observer would find that the flux from the star would 
vary with a frequency of mJl — cu, slightly less than the spin frequency of the star for m = 1. 



2.2.1. A Degnerate Layer 

If the mode resides in the degenerate ocean lying beneath the burning layers, the fre- 
quency of the mode is given by Eq. 15 of Bildsten et al. (1996) or Eq. 2 as a — 0, 

/A,n.D-- = 2.37Hz^^^^-j (^^j(l + 4.2n^) ^xj^ Kelvin modes 

(12) 

in the limit of rapid rotation, i.e. large q = 2Q/uj. The ocean is left after a Type-1 burst which 
burns the accreted hydrogen and helium directly to the iron group elements (Lewin et al. 
1993) and possibly beyond (Schatz et al. 2001); for simplicity A is set at 56 ; Bildsten et al. 
(1996) were concerned with the carbon ocean remaining after the steady nuclear burning 
associated with the more rapidly accreting "Z" sources. This frequency estimate assumes 
that the material in the excited layer is degenerate while the mode is excited. 

The temperature at the bottom of the burning layers at the peak of the burst is typically 
around 10^ K and drops to 10^ K as the burst subsides (e.g. Joss 1978). For the u — 1 and 
n — 1 westbound r— mode in the burning layer itself, the degenerate formula (Eq.l2 yields a 
frequency shift from 1.62 Hz at the peak temperature of 2.2 x 10^ K (Gumming & Bildsten 
2000) to 0.35 Hz at 10^ K. Even this modest frequency shift would require that the massive 
ocean be well coupled to the burning layers above. Because the mass of the ocean is much 
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greater than that of the burning layers, it is unhkely that the ingoing flux during the burst 
is sufficient to heat up the ocean to 10^ K. 

Modes with higher values of n will exhibit smaller frequency shifts — for n = 2 the shift 
is 45% smaller. While modes with m < — 1 will have larger shifts (but the variabihty will 
be found near |m| times the rotational frequency), the observed variation in the total flux is 
less. The Kelvin mode exhibits a frequency shift that is three times larger because A is nine 
times larger (see Eq. 5). The observed frequency of the mode decreases as the ocean cools. 



2.2.2. A Non-degnerate Layer 

The burning layers become sufficiently hot to lift the electron degeneracy (Gumming 
& Bildsten 2000), the frequency shift would be an order of magnitude larger. Additionally, 
because I have assumed that electron scattering dominates the opacity, the ratio of the gas 
pressure to the total pressure {(3) is constant through the non-degenerate layer where neither 
convection nor nuclear burning occur. If I take f3 equal to unity, and the temperature of the 
burning layer going from 10^ K to 10^ K while the mode is excited, the mode frequency will 
change by about eight Hertz much larger than is observed. 

However, a mode could naturally become trapped in the photosphere itself where the gas 
is typically much cooler. Lapidus et al. (1994) argued that sound waves can be trapped above 
the photosphere of the neutron star. Mdoes are trapped above the photosphere of the sun in 
the chromosphere between the temperature minimum and the steep temperature gradient at 
higher altitudes (e.g. Leibacher et al. 1982) It is natural to speculate that something similar 
could occur in the atmosphere of a Type-I burst. Lapidus et al. (1994) considered situations 
where (3 <^ 1 and the photosphere expanded to several stellar radii {i.e. radius expansion 
bursts). Here I will focus on gravity waves in the relatively thin photospheres of burst that 
do not experience radius expansion. 

Because /3 is constant through the photosphere and it is potentially observable, it is 
natural to use the condition of hydrostatic equilibrium to eliminate the temperature from 
the frequency estimate. 

. (1 - /,)^ . (1 - ffl jj^^^l^ (2.4 X lO'K)' -d 1 - ^ ^ £ (13) 

This substitution in Eq. 3 yields 

,(l-/9)'/'4-3/5 0.6V^%, , ,o 2X-1/2 
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where I have assumed that the photosphere spans an order of magnitude in density. Two im- 
portant features of Eq. 14 is that it depends on observable quantities and that the frequency 
drift diverges as the flux approaches the Eddington hmit. 

If the burst oscillation depicted in Fig. 2 of Strohmayer & Markwardt (1999) is taken 
as an example with the assumption that the burst peaked very close to the Eddington flux, 
the oscillation starts with F/F^dd ~ 0.9 and continues until F/Fsad ~ 0.1, which yields a 
frequency shift of 4.1 Hz (or 3.07 Hz including a gravitational redshift of 0.35) about twice 
the observed value. Strohmayer (1999) found oscillations whose frequency increased as the 
burst subsides. The excitation of a Kelvin mode could naturally explain this behavior. The 
spin-down portion of the oscillation begins when the flux is 0.6 of its peak value and continues 
until the fraction is 0.3 and its frequency increases by 1.4 Hz. The model predicts a slightly 
larger shift of 1.9 Hz (again with the gravitational redshift). 

It is crucial to emphasize that the estimates of the frequency shifts given in the previous 
paragraph are upper bounds for a layer of a given thickness in density. If the composition of 
the photosphere is not solar but helium, the frequency shifts would decrease by one-third. 
Less subtle is the assumption that the bursts reach the Eddington rate at their peak. For 
example if the two example bursts discussed earlier peaked at 80% of the Eddington rate the 
frequency shifts would be 1.43 Hz for the spin-up example - almost exactly the value found 
by Strohmayer & Markwardt (1999) - and 1.30 Hz for the spin-down case. In sources with 
radius expansion bursts, the peak burst flux can be calibrated and these predicted trends 
tested in detail to possibly give hints of the radii and gravitational redshifts of the underlying 
neutron stars. 

The observed burst oscillations have been found to have Q— values of several thousand. 
Since the frequency of the mode is boosted by a factor of several hundred, this is consistent 
with the mode having a Q— value on the order of ten. The observed value of Q is given by 
the product of the ratio of the observed frequency to the inherent frequency of the mode 
and the Q— value of the mode. Such moderate values of Q are not difficult to achieve in 
the neutron star ocean (Bildsten & Cutler 1995). The oscillations of interest suffer little 
dissipation on the timescale of the bursts. These details of the photospheric modes have not 
yet been determined. 
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3. Discussion 

The preceding sections have presented a model for the observed varabihty in Type- 
I bursts and its shift in frequency. Unhke the simple rotational modulation model, the 
excitation and evolution of r-modes in the neutron-star ocean can naturally account for the 
size and sign of the shift and the presence of variability even after the burst has enveloped 
the entire star. 

During the onset of the burst when only a portion of the star is hot, it is quite natural 
to account for the variability as a growing hotspot (Strohmayer et al. 1997b). Strohmayer 
et al. (1998) found that the pulsed fraction of 75±15% during the onset of the burst. The 
presence and growth of this hotspot sets up travelling modes in the ocean. Due to a match 
in timescales between the growth of the spot and the oscillation of the modes, modes with 
frequencies near a Hertz are preferentially excited. Modes with m = 1 result in the largest 
observed variability. From the observer's point of view, the flux varies at slightly less than 
the rotational frequency of the star. During the decay of the burst, the observed modulation 
is significantly lower ~ 15% (Strohmayer et al. 1998). The r— modes with the simple coupling 
considered in the letter can easily generate this amplitude of modulation. 

Although modes with higher (more negative) values of m may be excited, the resulting 
variability is much smaller. Since the burst begins in a particular spot and grows, one would 
expect modes with odd values of m to be excited preferentially. The resulting variability for 
the m = 3 mode is a factor of 30-100 smaller than the m — 1 mode; this is well below the 
upper limits quoted by Muno et al. (2002a) for several sources. 

Miller (1999) discovered that 4U 1636-536 exhibits oscillations at 580 Hz as well as much 
weaker oscillations 290 Hz; on the basis of other observations (QPOs) and a model of the 
accretion flow (Miller et al. 1998), Miller argues that the latter is the spin frequency of the 
star and the former is its first overtone. He associates the 580 Hz signal with the presence 
of two nearly antipodal hotspots burning on the surface (Weinberg et al. 2001). Recent 
observations of burst oscillations from SAX J1808. 4-3658 at the known spin frequency of 
about 400 Hz clarifies this issue (Chakrabarty et al. 2003). The difference in the frequencies 
of the QPOs is only about 200 Hz, so in other sources where the burst oscillation frequency 
is twice the QPO difference, the burst frequency may also be the spin frequency of the star, 
so the |m| = 1 mode is indeed the correct one to consider. 

On the other hand, if 290 Hz is indeed the spin frequency of 4U 1636-536, in these 
particular bursts, anomalously, the m = 2 mode is excited more strongly that the m = 1 
mode. If both modes have the same value of n (the number of radial nodes), a prediction of 
the model would be that the frequency shift would be twice as large for the 580 Hz variation 
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as for the 290 Hz variation. The simple rotational modulation model predicts a similar trend. 
However, the value of ui could be chosen by the evolution of the burning front, so the two 
modes could correspond to different values of n with the nearly same value of a;. In this 
case, both oscillations would vary in frequency by the same amount. 

This brings up the important question of why the burst only excites a mode with 
a particular number of radial nodes. If this were not the case, power would appear at 
several frequencies near the spin frequency of the star and these frequencies would all evolve 
according to Equation (3). Perhaps it is not surprising that only a mode with a n = 1 is 
excited since the excitation of higher modes requires additional shearing in the middle of the 
excited layer. The burst is associated with activity only in a thin layer at the top of ocean 
or at the bottom of the photosphere. 

This model makes several predictions. In principle eastbound modes may be excited. 
This would result in an observed frequency decrease as the atmosphere cools after the burst 
about three times larger than that caused by the westbound modes (e.g. Strohmayer 1999). 
More generally, several modes could be excited. If several modes were excited, one could 
in principle derive properties of the neutron star since the spacing of the modes is well 
understood. These higher order modes would typically exhibit less variability since the 
variability is proportional to q"^/^ or uj^/^. 

Furthermore, since the frequency of the mode changes by a factor of several as the 
surface layers cool, one would expect the variability to decrease slightly during this epoch 
as well. Additionally, more quickly rotating stars should exhibit slightly less variability and 
more importantly similar frequency shifts, i.e. the frequency shift is not proportional to Vl 
in contrast with standard, angular-momentum conservation model. Gumming & Bildsten 
(2000) noted that the frequency shift rather than the fractional frequency shift is similar 
from source to source. Also, because the mode may still be excited after the ocean cools to 
its equilibrium temperature, one need not expect the observed frequency of the oscillation 
to asymptotically approach the precisely same value from one burst to the next. 

Muno et al. (2002b) noted that the asymptotic frequency is stable to a few parts per 
thousand from burst to burst. For the photospheric model, the limit of the observed fre- 
quency as the burst flux goes to zero is indeed the spin rate of the underlying neutron star 
and should be constant from burst to burst. Muno et al. (2002b) also found that a small 
fraction of bursts exhibited multiple oscillation frequencies or spin-down episodes which this 
model exhibits. Most interestingly they found that the burst oscillations cease during radius 
expansion episodes. During a radius expansion episode the frequency of the photospheric 
modes diverges so the photospheric model argues that no oscillations should be found during 
radius expansion. 
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Cumming & Bildsten (2000) argued that even a small magnetic field could play a dy- 
namic role in the frequency shift of burst oscillations. In their layer model, the shearing 
fluid could dramatically stretch the magnetic-field lines which could provide a significant 
resistance to the shearing motion and limit the expected frequency shifts. In this model, the 
fluid itself does not move a significant distance relative to the stellar surface so the magnetic 
field does not get stretched significantly. In principle a sufficiently strong magnetic field 
could change the mode frequencies (e.g. Morsink & Rezania 2002, see Eq. 1) but a detailed 
discussion of the effects of weak magnetic fields on these models is beyond the scope of these 
article. 

This letter proposes a model for Type-I burst oscillations in which the burst excites 
waves in either the degenerate ocean, the burning layer itself or the photosphere. The 
frequency in the rotating frame of the waves in the ocean and the photosphere are typically 
several Hertz which provides a good match with the deflagation timescale. The frequency is 
the burning layers is a factor of ten larger. 

Even though the ocean model does yield mode frequencies near the observed values, its 
thermal inertia is large so it is difficult to understand how its properties could change so 
dramatically as the burst subsides. The photospheric model best accounts for the observed 
frequency shifts. The predicted frequency shift depends on the potentially observable value 
of the ratio of the outgoing flux to the Eddington value, and it is close to the observed values. 

However, further work is necessary. The frequency estimates presented here were derived 
using several simplify assumptions. First, it was assumed that the mode frequencies scale 
simply with the Brunt- Vaisala frequency as one goes from a degenerate to a non-degenerate 
to a radiation-pressure-dominated regime. Second, the adiabatic estimates of Bildsten & 
Cutler (1995) and Bildsten et al. (1996) were used in the photosphere where a more realistic 
isothermal approximation would be appropriate. Third, the magnetic field of the neutron star 
could also affect the frequency of modes in the tenuous photosphere (see Bildsten & Cutler 
1995, for a discussion). A more detailed treatment of the physics of the photosphere of a 
neutron star during a Type-I burst could address all of these concerns as well as determine the 
photospheric structure of a neutron star during a Type-I burst and would be very worthwhile 
given the approximate agreement between these gravity-wave models and the observations. 

In general these gravity-wave models present the possibility of a family of modes being 
excited during the burst whose observed frequencies may increase or decrease as the ocean 
cools. Modes with more radial, azimuthal and latitudinal nodes tend to exhibit less variabil- 
ity when averaged over the observed portion of the stellar surface - the mode which lacks 
latitudinal and radial nodes and has only one azimuthal mode exhibits a significantly larger 
variability than the other modes, independent of viewing angle. Finally, the magnitude of 
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the frequency shift should be independent of the spin rate of the star, in contrast to the 
standard picture of Type-I burst oscillations. 
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